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Shape invariance is a powerful solvability condition, that allows for complete knowledge of the 
energy spectrum, and eigenfunctions of a system. After a short introduction into the deformation 
quantization formalism, this paper explores the implications of the supersymmetric quantum 
mechanics and shape invariance techniques to the phase space formalism. We show that shape 
invariance induces a new set of relations between the Wigner functions of the system, that allows 
for their direct calculation, once we know one of them. The simple harmonic oscillator and the 
Morse potential are solved as examples. 
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INTRODUCTION 

The phase space quantum mechanics formalism (also 
known as deformation quantization) uses complex func- 
tions acting on the phase space instead of operators act- 
ing on a Hilbert space. These functions are endowed 
with a novel multiplication rule, the *-product, which is 
non-commutative, associative, and hermitian. The for- 
malism maps operators A of quantum mechanics to com- 
plex functions A(x,p), and vice versa. Historically the 
fourth way of doing quantum mechanics, the deforma- 
tion quantization is an equivalent way of doing quantum 
mechanics. As Dirac noted [H] 

Two points of view may be mathematically 
equivalent and you may think for that rea- 
son if you understand one of them you need 
not bother about the other and can neglect 
it. But it may be that one point of view may 
suggest a future development which another 
point does not suggest, and although in their 
present state the two points of view are equiv- 
alent they may lead to different possibilities 
for the future. Therefore, I think that we 
cannot afford to neglect any possible point of 
view for looking at Quantum Mechanics and 
in particular its relation to Classical Mechan- 
ics. 

The goal of this work is to give a short introduction to 
the phase space quantum formalism, and to explore some 
of its implications as seen through the prism of the super- 
symmetric quantum mechanics (SUSQM). We show that 
shape invariance induces a new set of relations between 
the Wigner functions of the system, that allows for their 
direct calculation, once we know one of them. The paper 
is organized in three parts. In the first part we briefly 
present the main concepts of the phase space quantum 
mechanics formalism. In the second part we introduce 
SUSYQM and emphasize the role of shape invariance 
(SI) in the solvability of the system. In the third part 
we analyze some of the implications of SUSYQM and SI 



in deformation quantization and derive a recursion for- 
mula between the Wigner functions of the system. We 
calculate the simple harmonic oscillator and the Morse 
potential as a concrete examples. 

QUANTUM MECHANICS ON PHASE SPACE 

In the following, we will use the "hat" notation to sym- 
bolize quantum operators and the "non-hat" notation for 
the corresponding functions on the phase space. We will 
also consider only the one-dimensional case, where the 
phase space is simply (x,p). For the general case, see 

The star product 

To any quantum operator A we associate a phase space 
function A(x,p) using the Weyl's transform W(A) |34j . 
which in the coordinate \x) basis, reads 

W(A) = A(x,p) = h J dye^y(x+ r ^\A\x- T ^) . 

(1) 

All integrals run from — oo to +00, unless specifically re- 
stricted. The complex function A(x,p) = W(A) is also 
known as the Weyl's symbol of A. Without other addi- 
tional constraints, the reciprocal map W is not unique. 
A sufficient condition for the unicity of W^ 1 is obtained 
by choosing an operator ordering [3D]. We choose Weyl's 
ordering, which prescribes symmetrical ordered polyno- 
mials in x and p. I.e., W (xp) = (xp + px)/2. Thus, we 
write the reciprocal map as 

A(x,p) ee W-\A{x,p)) (2) 

= tA^ / dudvdxdpA(x,p)e l{ P- p)u+l(£ - x)v . 
(2tt)^ J 

With this definition, W^(x) = x and W^ 1 (p) = p. 
More generally, one can readily check that W~ 1 ((ax + 
bp) n ) = (ax + bp) n . 
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Next, we define the phase space product between the 
Weyl's symbols. This product should correspond to the 
product of the quantum operators in the Hilbert space. 
Groenewold [3T] and Moyal [25] showed that using Weyl's 
ordering assumption, the product of functions on phase 
space can be written as 



A(x,p) * B(x,p) = A(x,p) e : 



] B(x,p), (3) 



where the arrow indicates the direction in which the 
derivative acts. Using ([I]) one can prove that indeed 
W{AB) = A{x,p) * B(x,p). 

The Groenewold-Moyal product is also known as the 
star-product (^-product) of A(x,p) and B(x,p). It can 
be expressed in several equivalent forms, such as 

A(x,p)*B(x,p) = e^( a * d p'- a » d *')A(x,p)B(x',p'), (4) 
calculated at (x',p') — (x,p); or, using the Bopp shifts 



A(x,p)*B(x,p) = A(x+ y d^, p- y d x \)B(x,p\h) 



= A(x,p)B(x- y dp, p+ y d x ). 

The *-product can also be represented as an integral, 
such as the Fourier representation (3] [33] 



A(x,p) * B(x,p) = 



dx\dp\dx 2 dp 2 



(6) 



(nh) 2 
A(x 1 ,p 1 )B(x 2 ,p 2 ) 

xe -T:[p( x l- x 2)+x{P2-Pl) + (x2Pl-X 1 p2)} 

or as the alternate integral 23^ 
A(x,p) * B(x,p) = 7—^5- / dx\dpidx 2 dp2 



(^) 2 j 

A(x + xi,p + pi)B(x + x 2 , p + P2) 



xe 



We conclude this definition by emphasizing the cyclic, 
trace-like properties [T2] of the *-product: 

J dxdp A(x,p) * B(x,p) = J dxdp A(x,p)B(x,p) 

= JdxdpB(x,p)*A(x,pl8) 

Thus, we have constructed a consistent way of moving 
back and forth from quantum operators -acting on a 
Hilbert space, to complex functions -acting on the phase 
space 

W(AB) = W{A)*W{B), (9) 
W~ X (A * B) = W~\A) W~\B). (10) 



This novel product is in agreement with the product 
rules for the operators in quantum mechanics. It is non- 
commutative 



associative 



A * B ^ B * A, 



A*(B*C) = (A*B)*C, 



and hermitian 



A * B = B * A. 



where the bar denotes the complex conjugation. 
For example, using Bopp's shifts we obtain: 

x*p = [x+ —d p p= xp+ — 



(11) 
(12) 

(13) 
(14) 



p*x = (p^y<^ x = P^-y. (I 5 ) 

Then, the ^-commutator, or the Moyal bracket of x and 
p is 



[x, p]* = x * p — p* x = ih, 



(16) 



which is consistent with the canonical commutation rela- 
tion [x, p] = ih. We see that the Moyal bracket provides 
a homomorphism with commutators of operators. 



Moyal bracket and the correspondence principle 

Here we will briefly show that it is the Moyal bracket 
which provides a consistent way to quantization, and not 
the Poisson bracket, as conjectured by Dirac. 

The ^-product and the Moyal bracket are h deforma- 
tions [3 [22 [21] of the usual commutative (point-wise) 
product of functions AB, and of the Poisson bracket 
{A , B} P = d x Ad p B - d p Ad x B. We can write 



A* B = AB + 0(H) , 

hA,B},={A,B} P + 0(h 2 ). 
in 



(17) 
(18) 



Thus, the name deformation quantization. Property ( 18 ) 
implies that 



lim ±r[A,B]. = {A,B} P , 
ft-j-o in 



(19) 



showing that in the classical limit one recovers the Pois- 
son bracket. 

The Moyal bracket obeys Jacobi's identity 

[[AB]„C], + [[C I A]»,B]» + [[B,C]„4=0, (20) 

as well as the Leibniz rule 

[A,B*C\* = [A,B], *C + B* [A,C], . (21) 
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Thus, one can endow the space of Weyl's symbols not 
only with a Lie algebra structure with respect to the 
Moyal bracket, but also with an inner derivative. An 
analogous statement holds true for the commutator al- 
gebra of quantum operators, but not for the Poisson- 
bracket algebra. Groenewold and van Hove showed 
[2"T1 I3T] that there is no invertible linear map from all 
functions of phase space A(x,p), B(x,p), . . . to Hcrmi- 
tian operators in Hilbert space A,B,... that will preserve 
the Poisson bracket structure, i.e. W~ 1 ({A, B} P ) = 
±[W- x (A),W- l {B)]. Using Poisson's bracket as a 
starting point for quantization works only in the spe- 



cial cases where the series (17) terminates after the first 



order in h. This is true for functions at most quadratic 
in p and x. As counterexample, Groenewold [2TJ showed 
that the identically zero expression in Poisson brackets 



{x 3 ,p 3 } p + 3{xp 2 ,x 2 p} p = 0. 
in Dirac's quantization heuristics, becomes 



(22) 



in 1 J in 



ip 2 + p 2 x x 2 p + px 2 



= -3 k 2 , (23) 



thus exhibiting a deficiency of order Ti . 

Through the Weyl's invertible correspondence map ([l}, 
the Hilbert space of Hermitian operators endowed with 
the operator commutator, has as counterpart the alge- 
bra of Weyl's symbols endowed with the Moyal's bracket. 
Up to an isomorphism, the Lie algebra generated by the 
Moyal's bracket is the unique associative one-parameter 
deformation of the Poisson bracket pi 15] fT3 l [TB I [20 l 132] . 
This uniqueness is extended (up to an isomorphism) to 
the ^-product. Thus, it is the Moyal's bracket that gives 
the correct correspondence principle for the quantization 
scheme [U HQ EH] ■ 



The Wigner function 

Now that we have a coherent way of modeling quantum 
operators and their products, we turn our attention to 
modeling the quantum state. A pure state is described by 
the ket vector \tp), or equivalently, by the corresponding 
density operator p = IV'XV'I- We will consider here only 
pure states, because the generalization to mixed states is 
straightforward. Let us denote by P(x,p) the normalized 
Weyl's symbol of the density operator p. We have 



P(x,p) 



1 
1 

27 



(24) 



dye 



-tpy/ 



P{X, P ) = 



2tt 



2 



tiy 



2 ' 



dye'^^{x-^-)^{x+^-) . (25) 



P(x,p) is the celebrated Wigner function [33], and will 
play a central role in the deformation quantization tech- 
nique. Here are several essential properties of the Wigner 
function [H [53] 

(i) P(x,p) is real. 

(ii) 



dpP(x,p) = \*p(x)\ 2 = (x\p\x) (26) 
dxP(x,p) = \TP(p)\ 2 = (p\p\p) (27) 
dxdpP(x,p) = 1. (28) 



/ 



(hi) If P^(x,p) and P,p(x,p) correspond to to the states 
ip(x) and 4>(x) respectively, then 



2irh dxdpP^(x lP )P t p(x,p) . (29) 



dxtp{x)4>{x) 



The last property has two interesting consequences. 
If ip(x) = 4>{x) then 



J dxd P P 2 {x lP ) = ^, (30) 
and, if we choose ip(x) orthogonal to 4>{x) , we get 
dxd P P^(x,p)P c f > (x,p) =0. (31) 



Equation (31) implies that P(x,p) cannot be ev- 



erywhere positive. Because it can also take nega- 
tive values, the Wigner function is also known as a 
pseudo-distribution. 

(iv) P(x,p) is the only pseudo-distribution for which 
each Galilei transformation corresponds to the 
same Galilei transformation of the quantum me- 
chanical wave function 36i. I.e., if ip(x) i— > ip(x + 
a), then P(x,p) t— > P{x + a,p), and if ip(x) t— > 
exp(ip'x/h) ip(x), then P(x,p) n- P{x,p — p'). 

(v) If ip( x ) l— ^ VK - x )i then P(x,p) <-> P(—x, —p), and 
if Tp{x) i — y ip(x) then P{x,p) P(x, —p). 

To describe the time evolution of the system in the 



phase space, let us take the time derivative of (24), and 



use Schrodinger's equation together with its conjugate. 
We get 



Wt P ^ P) = 2^ 

= \lf 
in 

or 



dye 



-ipy i 



^\{Hp-pH)\*- n i) 



H. 



2-Kh 



), 



ih-P(x,p) = [H(x,p),P(x,p)]< 



(32) 
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where H(x,p) is the Weyl's symbol of the Hamiltonian 



H. Equation (32) is mirroring the time evolution of the 
density operator ih-j^p = [H,p\. The time evolution 
of the Wigner function can also be symbolically written 
using the sine notation @], which emphasizes the non- 
linear deformation involved by the Moyal bracket 



loss of generality we take 2m = 1,cj = 2. Then the 
Hamiltonian operator of the harmonic oscillator reads 



H =f + x 2 . 



(39) 



d_ 
di 



P(x,p) 



d d d d 
dx dpi dxi dp 



having the corresponding Weyl symbol 

H[x,p) =p 2 + x 2 . (40) 
H{x,p) P(xi,p^ e ^ ug wr j^ e ^-eigenvalue equation (34 1 for this case. 



calculated at (xi,pi) = (x,p). 

For stationary states ^P(x,p) — 0, hence the Hamil- 
tonian ^-commutes with the Wigner function 



[H(x,p),P(x,p)]. = 0. 



(33) 



In the stationary case one can get more constraints on 
the Wigner function. It we apply the Weyl's transform 
on Hp — Ep, we obtain the *-eignevalue equation [TT1IT5] 



H(x,p)*P(x,p) = EP(x,p), 



(34) 



where E is the energy of the system. Hermiticiy implies 
that pH = Ep, hence the symmetrical relation holds true 
as well 



P(x,p) * H(x,p) = EP(x,p) . 



(35) 



Note that if Pe and Pe< correspond to the eigenenergies 
E and E' , then due to the associativity of the *-product, 
we have 



P E * H * P E > = E P E * Pfj = E' P E * Pe' 



(36) 



If E 7^ E' , it follows that Pe*Pe> = 0. In the general case 
we have Pe * Pe 1 = \pE Pe') I \2ttK) 2 which yields the 
orthogonality-idempotence relation of Wigner functions 



Pe * Pe- 



°E,E' 

2itTi 



Pr 



(37) 



In deformation quantization, P(x,p) plays an analogous 
role to the probability density function in classical sta- 
tistical mechanics. Namely, the average of A in state l^) 
is given by [35) 

(ip\A\ip) = (A(x,p)) = / dxdpP{x,p) * A(x,p) . (38) 



Note the resemblance with statistical mechanics. How- 
ever there is a major difference in the case of the de- 
formation quantization formalism: the function P(x,p) 
is not a probability distribution in the statistical sense, 
because it can take negative values. 



The harmonic oscillator 

To exemplify the concepts introduced so far, let us con- 
sider the case of the simple harmonic oscillator. Without 



We obtain 

(p 2 + x 2 )*P(x lP ) = EP(x,p), (41) 
which, using Bopp's shifts representation ^ becomes 



V - -wOx 



P(x,p) = EP(x,p). 

(42) 



Performing the algebra we arrive to 



P(x,p) = 0. 



P 2 +x 2 + ih(xd p - pd x ) -^(d; + d 2 ) - E 

After separating the imaginary and the real parts, we get 
{xd p -pd x )P(x,p)=0, (43) 



p 2 + x 2 - h -{dl + d 2 x )-E 



P(x,p) =0. (44) 



The symmetry of equation (43 1 indicates that P(x,p) 
depends effectively of only one variable, which is a sym- 
metric combination of x and p. Thus, in ( |44[ ) we make 
the change of variables t = 2(p 2 + x 2 )/T% and write 
P(t) = e- l l 2 L{i). We get 



td 2 + (l-t)d t + 



(-- 




\2h 





L(t) = . 



(45) 



For L(t) to be normalizable, the zero derivative term of 
( 45 1 must be a positive integer 



E 
2ft 



72 = 0,1,2, . 



(46) 



This assures that the series solution of L(t) terminates at 
a given rank, hence L(t) is finite and e - '/ 2 L(t) is normal- 



izable. But condition (46) yields exactly the quantization 



formula for the energy of the harmonic oscillator (oj = 2) 

r 



E„ 



2ft 72 



72 = 0,1,2,. 



(47) 



Consequently, equation (45) subject to the constraint 



( 46 1 , becomes the differential equation of Laguerre poly- 
nomials L n (t) 9J. Thus, we obtain the analytic expres- 
sion for the Wigner functions P n (x,p) of the harmonic 
oscillator 



(-1)" p 2 +* 2 
P n (x,p) = - —e * L rr 



p 2 + X 2 

ft/2 



, 72 = 0,1,2,... 

(48) 
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Figure 1: The Wigner function Ps(x,p) for the simple har- 
monic oscillator. 



where L n is the n-th Laguerre polynomial. The harmonic 
oscillator is one of the few systems where the Wigner 
functions are completely known analytically pQ . In figure 
([!]) we illustrate the Wigner function for n = 5, m = 1/2, 
and u) = 2. Note the circular x ■H- p symmetry of the 



solution, as reflected by the equation ( 43 ) 



Finally, let's analyze the algebraic solution induced by 
the factorization of the Hamiltonian in terms of creation 
and annihilation operators a + and a. This method will 
naturally segue into the factorization of a general Hamil- 
tonian, which is the crux of the SUSYQM techniques. 
We write the Hamiltonian as 



H = 2h(a + a+1/2) 



(49) 



where a — (ip+x)/V2H and a + = (a)t = (—ip+x)/V2~h. 
Then [a, a + ] = 1, and after algebraic manipulations we 
obtain 

H\n) = 2h(n + 1/2) \n) , n = 0,1,2,... (50) 
a+ |n) = y/n + 1 \n + l) (51) 
a\n) = y/ri\n-l). (52) 

The harmonic oscillator wave functions are obtained by 
projecting the recursion relations 



Initio), 



(53) 



onto the \x) basis, where the starting point is given by 
o|0) = 0. 

In the phase space picture these relations become 

H(x,p) =p 2 +x 2 = 2Ti{a + * a + 1/2), (54) 

where the annihilation and creation functions a(x,p) and 
a + (x,p) are 



A simple calculation shows that their Moyal bracket gives 

[a,a+], = 1. (56) 

Then, the *-eignevalue problem becomes 

H(x,p)*P n (x,p)=2h(n + 1/2) P n (x,p), (57) 

where the Wigner function P n (x,p) is determined from 
the recursion relations 12J 



P n {x,p) = - {a+*) n P (x,p) (*ay 
nl 



(58) 



using as starting point 

a*P {x,p) = 0. (59) 
A direct calculation 12: shows that the algebraic method 



recovers the known form for the Wigner function ( 48 1 



This algebraic technique will be extended to a gen- 
eral Hamiltonian, by applying SUSYQM concepts to the 
phase space formalism. 



SUPERSYMMETRIC QUANTUM MECHANICS 

In this section we give a brief introduction to super- 
symmetric quantum mechanics (SUSYQM) techniques 
which we will use mostly as reference for the next sec- 
tion. For more details see PS EH QB [37]. SUSYQM 
provides an elegant and useful prescription for obtaining 
closed analytic expressions both for the energy eigenval- 
ues and eigenf unctions for a large class of one dimensional 
potentials. 

Given the Hamiltonian H = f— + V(x) with the 
eigenenergies E n , let us shift its energy spectrum by Eq, 
such that the new ground state energy becomes zero, and 
denote the new Hamiltonian by , and its eigenenergies 
hyE~: 



H-=^ m +V_(x); #_|n; -)=££>;->, £ "=0. 

(60) 

H- can be factorized I2U as 



where 



ip 

f2^L 



H- = A+A , 



+W(x), i+ = (i)t = 



-ip 



'2m 



(61) 



+ W(x). (62) 



The function W is known as the superpotential and is 
related to V-(x) by 



a = (ip+x)/ V2h , a+ = (a) = (-ip+x)/V2h . (55) 



V-(x) = W 2 (x) 



'2m 



W'{x) 



(63) 
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Equation (63) becomes a Riccati equation when pro- 



jected on the \x) basis, thus allowing for finding W(x) 
A sufficient condition for E$ = is given by 



A 10;-) =0. 



(64) 



which yields the explicit action of W(x) on the ground 
state eigenvector | 0; — ) 



W(x)\0;-) = - 



ip 

1 2m 



10;-) 



(65) 



If we project (65 1 on the \x) basis, we can express 1U(;e) in 



terms of the ground state eigenfunction ip (x) — (x\0; 

as \m 



W(x) = -- 



ft A '(x) 

/ 2m ipo(x) 



(66) 



By interchanging the operators A + and A we generate 
a new Hamiltonian H + = A A + , which corresponds to a 
new potential V+(x) 



H + = 



P 

2m 



+ V+(x); V+(x) = W 2 (x) 



'2m 



W'(x) , 



(67) 

Let us denote by E+ its eigenenergies: H + \n;+) = 
E+ \n; +). The two Hamiltonians if_ and H + are known 
as super symmetric partner Hamiltonians. Their eigen- 
vectors and eigenenergies are related 



E n - E n+l ! E - > 



:A\n + l:-) 



E 



n+l 



|n + l;-) 



A+\n;+). 



(68) 
(69) 

(70) 



Thus, with the exception of the ground state, the su- 
persymmetric partner potentials H— and H + share the 
same energy spectrum (isospectrality), and have inter- 



connected eigenvector sets. Formulas (68 - 70) together 



with (64) represent a quick algebraic way of completely 



finding the spectrum and eigenvectors of H + if we know 
the spectrum and eigenvectors of H- . In this way, if H- 
corresponds to a very simple simple system, we have an 
elegant solution for solving more complicated cases asso- 
ciated to H + , 

There is a special class of partner Hamiltonians for 
which potentials V-(x) and V+(x) obey an additional 
constraint known as shape invariance. Shape invariance 
states that partner potentials are similar in shape 



V + (x;a ) = V_(x;ai) + R(a ) . 



(71) 



differing only by a set of parameters (modeling the 
strength of the interaction) ag, ai = /(^o), and a con- 
stant term i?(ao), independent of x. Shape invariance is 



a solvability condition that allows in principle to com- 
pletely solve the system, as it will be shown shortly. For 
this, let us rewrite the shape invariance condition as 

H + (x,p;a ) +g{a Q ) = H-(x,p;ax) + g{a x ) , (72) 

where R(a ) = g(ai)—g(a ). Note that we do not have to 
know the spectrum of one partner Hamiltonian to know 
the other. Since Hamiltonians in (72) differ by a con- 



stant, their eigenvalues differ by the same constant, and, 
up to a global phase factor, both have the same set of 
eigenvectors. Therefore, for all values of n we have 

E+(a Q ) = E~(ai)+g(ai)-g(a ) ) (73) 
\n,a ;+) = |n,en;-). (74) 

For normalizable ground states (unbroken SUSY), the 
ground state energy of H_(x,p;ai) is zero, Eq (ai) = 0, 
for each iteration of the parameter a, = /(aj_i). By suc- 
cessively using the shape invariance (|73[) and the isospec- 



trality (68) of the partner Hamiltonians, in conjunction 



with the unbroken supersymmetry condition, we get 



E o fao) 
E i(a ) 



= 0. 



= £^(ao) (isospectrality) 

= Eq(o,i) + g(ai) — g(do) (shape invariance) 

= g(ai) ~ g(a ) (unbroken supersymmetry) 



K ( a o) = g(a n ) - g(ao) ■ 
Similarly, 



n,a ; 



i+(a o )---i + K-i)|0,a„;-) 



1/2 



(75) 



(76) 



All eigenenergies and eigenvectors of the Hamiltonian 
can be determined iteratively by this algorithm. 
Thus, SUSYQM and shape invariance determine the en- 
tire spectrum of the system without any need to solve 
complicated differential equations. For more details see 
|18) and the references therein. 

We are now ready to apply these ideas to the phase 
space formalism. 



SUSYQM AND SHAPE INVARIANCE IN PHASE 
SPACE 



SUSYQM on the phase space is obtained by apply- 
ing Weyl's map to the operatorial framework discussed 
above. Thus, the Weyl symbols of the partner supersym- 
metric Hamiltonians and H+ are 



H-(x, P ) = ?- + V-(x) = A + (x,p)*A(x,p),(77) 
Zm 



H+(x,p) 



2m 



+ V+(x) = A(x,p) * A+(x,p) , (78) 
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where the annihilation and creation functions are given 
by 



A(x,p) = ^+W(x); 
2m 

A+(x,p) = A(x,p) = -^+W(x). 

2m 

The supersymmetric partner potentials are simply 

h 



V t (x) = W 2 (x)t 



'2m 



W'{x) . 



(79) 
(80) 

(81) 



The last formula can be obtained either by Weyl's map- 
ping the quantum operators V^(x), or by direct cal- 



culations from ( 77 ) and ( 78 1 using the definitions ( 79 



80), and the property of Moyal's bracket [W^(a:),p]* = 
ihd x W{x). 

The ^-eigenvalue problem for the supersymmetric part- 
ner potentials H_(x,p) and H+(x,p) written in terms of 
their corresponding Wigner functions reads 



H-(x,p) * P n (x,p) = E n P n (x,p) 
H+ (x,p)* P+ (x, p) = E+ P+ (x, p) . 



(82) 
(83) 



where 



2tt 



dye 



-ipyl 



2 ' 



(84) 

are the n-th excited state Wigner functions of the super- 
symmetric partner Hamiltonians H^(x 7 p). The neces- 
sary condition for the ground state energy Eq of iJ_ (x, p) 
to be zero, is 



A(x,p) *P (x,p) = 0. 
This yields the ^-product equation 



W(x)*Pv(x,p) 



P* P Q (x,p) , 



V2m 

which in terms of Bopp's shifts becomes 

—i 



(85) 



(86) 



W[x+~d p ) P ~(x,p) 



2niK P--d x }P»{x,p). 

(87) 



By applying the Weyl map on equations ( 69 ) and 



(701 we obtain [TT] the connections between the *- 



eigenfunctions of the supersymmetric partner potentials 
H_(x,p) and H + (x,p) 

P n(x,p) = -J—A(x,p) * P~ +1 (x,p) * A + (x,p) 



E 



n+1 



E n 



A+(x,p)*P+(x,p)*A(x,p). (89) 



The only constraint we imposed so far was that 
the Hamiltonians H^(x,p) are supersymmetric part- 
ners. This leads, besides the isospectrality, to the *- 
product connection between the corresponding Wigner 



functions P n f (x,p). It is almost the "carbon copy " of 
the SUSYQM case. Next, we explore the implications 



of the shape invariance. In phase space, equation (72) 
becomes 



H + (x,p;a Q ) +g(a Q ) = H^(x,p;ax) + g{a 1 ) 



(90) 



where ao and a± = /(ao). Because the Hamiltonians 
in (90) differ by a constant, they have the same set of 
*-eigenfunctions. Hence, up to a sign factor (because 
Wigner functions are real), shape invariance implies 

Pn( x >P\ a o) = Pn( x ,P', a i) > (91) 
and the energy spectrum 

E n( a o) = g(a n ) - g(ao) , with a n = /"(ao)- (92) 

In addition, the shape invariance together with the su- 
persymmetric condition, leads to a new recursion formula 
among the Wigner functions: 

P n (£,p;a ) = A^(x,p;a ) * ^_ , _ , * A{x,p;a ) 



The proof is immediate: 



P n ( x ,Pl a a) = A+ (x,p;a )* 



E n (a ) 



-P+_i(x,p; op) 



(93) 



A(x,p; a ) 



= A^{x,p:,a )* — — * A{x,p;a ) , 

E n (ao) 

where the first equality follows from SUSYQM, and the 
second one from shape invariance and isospectrality. We 



can now iterate ( 93 ) to obtain the expression for a gen- 



eral P n (cc,p;ao) starting from the ground state Wigner 
function Pq(x,p; a„). We have 



P n {x,p]a ) = A + (x,p;a ) * ■ ■ ■ * A+(x,p;a„_i) 
Po(x,P;a n ) 

*■■■* A(x,p; a ) 



A(x,p;a n 



(94) 



This recursion formula together with ( 85 ) as its starting 



point, determines, the entire set of Wigner functions for 
shape invariant systems. However, for concrete exam- 
ples, the multiple *-products can become prohibitively 
difficult to calculate. The harmonic oscillator is one a 
the few cases where we can determine analytically the 
entire sequence of Wigner functions. 



The Morse potential 

As an example, let us consider the one dimensional 
Morse potential. The Morse oscillator is a good ap- 
proximation of the oscillatory motion of the bi-atomic 
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molecules [5] [57] • Without loss of generality, we take 
h = 1 and 2m = 1. Then, the Hamiltonian reads 



H- (x, p: a) = p 2 + V- (x, a) 
where the Morse potential 



(95) 



V-(x,a) = a 2 +b 2 e 



2-2sx 



26(a+l/2)e" 



a,b,s > , 
(96) 

can be generated as V- — W 2 — W from the superpo- 
tential [10] 



W(x; a) = a — be 



(97) 



The corresponding annihilation and creation phase space 
functions are 



A(x, p) = ip + a — be 



A + (x, p) = —ip + a — be 



(98) 

As we know [TU], Morse potential is shape invariant. The 
relevant parameters of the model are ag = a and cii = 
/(do) = a — s respectively. Thus 



H + (x,p;a)+g(a) = H_(x,p;a-s)+g(a-s) ; g(a) = 



and consequently, the energy spectrum is 
E~(a) = a 2 — (a — ns) 2 . 



(99) 



(100) 



Let's find the ground state Wigner function P (x,p) us- 
ing 



A(x,p) *P (x,p) = 0. 



(101) 



Writing the *-products in (ToTJ as Bopp's shifts, we ob- 
tain 



i p 



l 7? 
2 



a — be 



P o -(x,p) = 0, 



which becomes a mixed differential/finite-difference 
equation 

{a + ip)P^{x,p) = be- sx Po(x,p-is/2) 

- ±d x Po(x,p). (103) 

The above equation mixes the partial derivative with 
respect to one variable with a shift with respect to 
the other variable. This hints to a well-known [26] 
differential/finite-difference formula for the modified 
Bessel functions K u (x) 



2d x K v {x)+K v _ 1 {x)+K 1+v {x) = Q. 



(104) 



The exponential term in the r.h.s. of equation ( 103 ) sug- 



gests that K„(x) depends on x via an exponential term. 
The shift in momentum variable suggests that p appears 



in the index v of K v (x). Therefore we substitute in ( 103 ) 
the following ansatz 



Po(x,p) = Cie~ C2X K aip) (C 3 e~ sx ) 



(105) 



where the constants C\, C2, C3, and the function a{p) are 



to be determined. Plugging (105) into (103), after some 



algebra, we find that ( 103 1 is identically satisfied for 

2ip 



C2 = 2a , C3 = — ; and a(p) 



(106) 



for any value of Ci, which is constrained by the nor- 



malization of the Wigner function (28). We obtain 



C\ = ^(2b/a) 2a / s . Thus, the ground state Wigner func- 
tion for the Morse oscillator is 



P (x,p;a) 



2b 



its \ a 



2a/ s 



e~ 2ax K. 



2b 

2i P /s \ — e 



. (107) 

Finding the next excited states Wigner functions is 



now straightforward using the shape invariance ( 93 ) . For 
example, let us find Pi(x,p;a) 



„ , x . n Pn (x, p: a — s) 

Pf (x,p; a) = A+ (x, p; a) * ^ z 1 

E 1 (oo) 



A(x,p; a) 



[a-ip-b e- sx ] * \ {x f a $ *[a + ip-b e~ ax ] 
a 1 — (a — s) z L 



1 



a-i[p--d x j - be- sl - x+ ^ a * ] 
a + i(p+^j -be-< x -^ 



After some algebra, we arrive to 



Pp {x,P\a~ s) 
2as — s 2 



P l {x,p;a) = a{x)K v {y)- p{x)[K v -i{y) + K v+ i(y)\ 

(108) 



where & 



— a — 2ax 



(102) a(x) 



2^e 



< and 

2+ 2a. 



\a-s) 



2n(2a - s)s 2 

^a/s^ e (~2a+s)x I __b_ 



[4b 2 + e 2sx {2a- s) 2 



P(x) 



In figure ([2| we illustrate the Wigner function P 1 (x,p) 
for the Morse oscillator. The graph is cut on the vertical 
axis to emphasize the details near the origin. 



CONCLUSIONS 

We have shown that quantum mechanics can be ex- 
pressed equivalently in the language of "normal" func- 
tions defined on the phase space, endowed with the *- 
product. The Moyal bracket corresponds to the operator 
commutator, and one can use it as an heuristic tool for 
quantization. In the classical limit (h — > 0) the Moyal 
bracket becomes simply the Poisson bracket, and the *- 
product, becomes the normal commutative product of 
functions. In this picture, quantum mechanics appears 
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Figure 2: The Wigner function P\{x,p) for the Morse oscilla- 
tor, for a = 5, b = 1, and s = 1. 



as a deformation of the classical mechanics, with the de- 
formation parameter h. Supcrsymmetric quantum me- 
chanics induces additional relations between the Wigner 
functions, while shape invariance exposes a simple ap- 
proach to recursively obtain the Wigner function of a 
system, starting from its ground state Wigner function, 
as exemplified for the Morse potential. 
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